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We revise a previous version of this paper [2| proposing a method to transfer emittances from the 
transverse planes to the longitudinaL It turned out the system showed non-Hamiltonian behavior 
because the longitudinal component of the transverse electric field at ends of the transverse resonator 
were not included in the model. Including those fields the 6D emittance is conserved. We review the 
mechanism behind the emittance conservation. It is still possible to produce an emittance transfer 
with the RF quadrupole field setting the tunes on synchro-betatron resonance for a short period of 
time then reducing the transverse RF field to zero. We found some coupling between the transverse 
and the longitudinal planes is always present, even there in no special device in the accelerator to 
provoke it. The role of the transverse RF field is to modify the transverse tunes quickly putting the 
system to synchro-betatron resonance for a short time. The emittance tends to flow from the plane 
with bigger emittance to the smaller one. 



I. INTRODUCTION 

In accelerators one of the most important parameter of 
a particle beam is its emittance, in many cases reducing 
the emittance is highly important. The common method 
to achieve small emittance is beam cooling. Cooling 
methods reduce the volume of the 6 dimensional phase- 
space. There are several types of cooling methods, each 
having its limitations, and in many cases non of those 
method can be used, or they are not practical. In certain 
cases it is more important to have a small emittance in 
one or two planes of motion than the others. Currently 
known methods , Q for complete emittance exchange 
between subspaces are subjects to certain restrictions. In 
accelerators where the beam transfer is symplectic and 
linear there are constraints which prevents the contin- 
uous transfer of emittances between different planes of 
motion, only a complete exchange is possible. This was 
first proved by E.D Courant M and it is known as the 
emittance exchange theorem An other method @ 
uses coupling resonances in a special ring to change the 
projections to subspaces. Methods using coupling reso- 
nances can produce partial emittance transfers between 
the planes of motion. In the previous version of this paper 
[H we presented a new method believed to produce non- 
Hamiltonian effects. Unfortunately the non-Hamiltonian 
behavior was a result of a missing longitudinal compo- 
nent of the transverse electric field at the ends of the res- 
onator which wasn't included in the model. We re-cite 
the idea described in the previous version I'j and discuss 
it with the correct model including the end fields. From 
time to time attempts appear to reduce the 6D emittance 
by external EM fields. K.W.Robinson has proved in the 
appendix of his paper 0] , that the synchrotron radiation 
damping partition numbers can not be changed by a gen- 
eral EM field. One aim of this revision is to recall this to 
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who might try similar methods. Another reason is that 
it is possible to produce a transverse-longitudinal emit- 
tance transfer with the device described below, although 
this emittance transfer is subject to all restrictions im- 
posed by symplectic conditions. An important difference 
compared to the results obtained in the previous version 
of this paper [l| is that the emittance tends to flow from 
the plane having the larger emittance to the plane with 
smaller emittance. Existing proposals 3, fo- 

cus on application for free electron lasers. Our method 
doesn't require a separate ring for the emittance transfer, 
the transverse resonator can be placed in an existing ring. 
Its role is not to introduce a coupling between the longi- 
tudinal and transverse planes, but to bring the transverse 
tunes to synchro-betatron resonance conditions and mov- 
ing away quickly when the required amount of transfer 
has happened. The simulations showed there is always 
some coupling between the longitudinal and transverse 
planes, even when the transverse resonator is not present 
and the accelerating cavity is placed in a dispersion free 
region. 

II. THE BASIC IDEA REVISED 

We re-cite the idea presented in the previous version of 
this paper [l| and also the reason why it was believed to 
produce non-Hamiltonian effects. Then we describe the 
mechanism preventing us from using the transverse res- 
onator that way. Let's consider a particle circulating in a 
storage ring, that contains at least one accelerating cav- 
ity with some RF voltage, and another RF resonator that 
generates a transverse quadrupole electric filed. This sec- 
ond electric field can be given as 

'E.{x,y,t)=Eg{ix~]y)cos{ujt + (j)) , (1) 

where Eg — dE/dx is the gradient of the electric field 
m [V/m^] and i, i are the Cartesian unit vectors. The 
electric field in eq. (U) can be derived from the vec- 
tor potential, the magnetic field obtained from the same 
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vector potential is zero. 
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(3) 

(4) 



The equations above give a good model of the fields in- 
side the transverse resonator (TR from here). There are 
however longitudinal components of the electric field at 
the ends of the TR not included in the model of our first 
attempt. These end fields can be modeled as a longitu- 
dinal kick giving an energy change to the particle 
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It is straightforward to obtain eq.([n]). If we choose the 
potential to be zero outside the TR, a particle traveling 
through the end field will get a longitudinal kick corre- 
sponding the the potential difference, which is propor- 
tional to the square of its transverse position. The direc- 
tions of the end fields are opposite at the entry and the 
exit. 

Lets consider first the effect of the transverse electric 
field inside the TR. The direction of the field points in- 
ward in the horizontal plane and points outward in the 
vertical plane. Fig. [1] and Fig. [2] might help visualize 
this field. Imagine a beam diverging in the horizontal 
plane and converging in the vertical. A positively charged 
particle traveling in such a beam looses kinetic energy 
through the RF electric field, because its motion most of 
the time has transverse components which are opposite to 
the direction of the force acting on it. There is an energy 
flow outward from the system each time a bunch of par- 
ticles goes through the TR. The energy loss due to the 
transverse electric field is compensated to some extent 
by the energy change due to the end fields. At the elec- 
trostatic limit the energy changes due to the transverse 
and longitudinal components cancel each other exactly. 
There is a frequency and length of the TR at which a 
the phase of the RF field changes tt radians while a syn- 
chronous particle travels through the TR. In this case 
the energy change is purely due to the transverse elec- 
tric field. The amount of the energy loss depends on the 
angle of the particle's trajectory and also on its longi- 
tudinal position. In the accelerating cavity the particle 
gains or loses energy depending only on its phase with 
respect to the synchronous particle. This later energy 
change is strictly a function of phase, while the energy 
change in the TR is a function of the transverse trajec- 
tory inside that field and also the phase of the particle 
with respect to the transverse RF field. It is not easy 
to see why the energy fiowing outward the system at the 
TR equals exactly the energy change in the accelerating 
cavity averaged over some period of time. 

The Panofsky-Wenzel theorem already indicated 
why an EM field can not reduce the 6D emittance. Later 
Hereward ^ found that for any electromagnetic field in a 



cavity the following equation must hold as a consequence 
of Maxwell's equations. 



dz dx 



(6) 



where z is the longitudinal position deviation from the 
synchronous particle. Eq. ([6]) shows that the force acting 
on the particle is a conservative one. In the appendix of 
[3] Robinson obtained the same equation and also gave a 
detailed derivation of why a general EM field can't change 
the synchrotron radiation damping partitions. The same 
mechanism is valid for our case. We re-cite here only the 
main points from The equation of the longitudinal 
motion ignoring the synchrotron radiation can be written 
as 



dE 

dN 
dz 
dN 
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Where e is the energy change per turn due to the EM 
fields in the ring, I is the path length, E is the energy of 
the particle and z is the longitudinal position deviation 
from the synchronous particle. There are two terms in 
the equations above which can produce damping and it 
is shown in Q that these two terms cancel each other as 
a consequence of eq.® 
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dz 
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In other words, the longitudinal damping produced by 
the energy loss depending on the energy is compensated 
by the anti damping due to the path length dependence 
on longitudinal position or vice versa. The path length 
depends on the longitudinal position because it is pro- 
portional to the transverse deflection introduced by the 
TR field, which is itself dependent on the longitudinal 
position of the particle. 

An other similar case has been analyzed in , where 
a dipole mode cavity was placed in a dispersive region. 
They came to the same conclusion, energy loss depen- 
dending on the energy is compensated by the the path 
length change introduced by the dipole mode cavity. 

In both case using the TR or a dipole mode cavity, 
the force introduced by the EM field is a conservative, 
therefore Liouville's theorem applies. The 6D emittance 
stays constant. Even with this constraint it is possible to 
introduce correlations between the planes of motion and 
change the projected emittances using synchro-betatron 
resonances. 



III. TRANSVERSE ELECTRIC FIELD 
GENERATION 

It isn't an aim of this paper to present a detailed design 
of the apparatus, but to demonstrate that the method is 
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FIG. 1: One possible design of a TR to generate quadrupole 
electric field 
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FIG. 2: The electric field between the the rods. 



IV. SIMULATIONS 



The method 



A Java code has been developed from scratch for 
the simulations. It was believed the TR produce non- 
Hamiltonian effects, and that was the main reason to 
develop the code and not using standard tools. The ef- 
fect we are interested in accumulates over thousands of 
turns, so it was important to avoid any approximations 
and have a full understanding of what the code does. 
Classical methods like Runge-Kutta-Cash-Karp with au- 
tomatic step size adjustment was found to be precise 
enough for the number of turns simulated. The source 
code is freely available from the author on request. Five 
type of elements have been implemented: Sector bending, 
magnetic quadrupole, RF accelerating cavity, drift space 
and the TR. Simulations used some combination of these 
elements. All elements were treated with the hard-edge 
model, expect the TR where the end fields turned out to 
be essential. The numerical precision has been verified 
by performing time reversed tracking. The initial condi- 
tions of the forward tracking has been re-obtained with 
6-7 digit precision after time reversed tracking (20000 
traversals of the FODO cell). 



1. Sector bending magnet 

The trajectory in the sector bending magnet is treated 
in a geometric way in a global coordinate system, where 
the x,y plane is the bending plane and the z coordinate 
is the vertical one. The particle arriving in the bending 
magnet moves on a circle with a radius determined by 



P = 



Pxy 



(9) 



feasible and the electric field of the TR can be practi- 
cally realized. The electric field is similar to the field 
of an RFQ, but the rods have no undulations. The fre- 
quency of the TR must be relatively low, because during 
the time a bunch travels through the TR the RF voltage 
should not change sign. One possible structure is de- 
scribed in [llj . This type of design has the advantage of 
having a resonant frequency inversely proportional to its 
length. The RFQ in [ll[ without vane modulation has 
been modeled with Microwave Studio [l^. The struc- 
ture is shown on Fig. [TJ The diameter of the cylinder is 
1 m, it's length is 2.5 m, the diameter of the rods are 
8.72 cm. The electric field between the rods has the best 
linearity if the ratio between the radius of the rods and 
their distance from the center is 1.14511. With these 
dimensions the structure resonates at 13.56 MHz. The 
frequency can be lowered making the structure longer. 
Fig. [5] shows the electric field from a view point on the 
beam axis. The magnetic field is mainly located outside 
the rods and has little effect on the beam. 



where 



Pxy = jmo^v^+u. 



(10) 



The center of the circle lies on the line perpendicular to 
the particle's trajectory at the entry. Then the inter- 
section of the trajectory circle and the exit edge is calcu- 
lated, that gives the x and y coordinates at the exit. The 
Vx ,Vy velocity components are rotated in a global coordi- 
nate system by the angle extent of the arc between the 
entry and the exit points, Vz is unchanged . The z coordi- 
nate at the exit is calculated as z = + VzL/vxy, where 
L is arc length between the entry and the exit points. 



2. Quadrupole 

The global coordinates and velocities are transformed 
to a local one, then the trajectory is calculated by solving 
the equations of motion numerically, and finally the local 
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coordinates and velocities are transformed back to the 
global coordinate system. In the local coordinate system 
X stands for the horizontal position y for the vertical and 
z for the longitudinal. The equation of motion is obtained 
as 



dp 

H 

B 



q{Bx v) , 

{Bgxit),Bgyit),0} , 
Bx 
dx ' 



(11) 



which gives 



z"{t) 



x"{t) + K^xz'{t) = 
y"{t)~Krnyz'{t)^0 
{xx'{t) - yy'{t)) = 



with 



(12) 



(13) 



Please note that there is a longitudinal component of 
the force, which is missing in the usual approximative 
formalism. The independent variable is time. The precise 
time when the particle exits the quadrupole isn't known 
in advance, it is found by successive approximations. 



3. Cavity 

The longitudinal component of the momentum is in- 
creased in a single step at the middle of the cavity's 
physical length, by the amount which corresponds to the 
energy change: AE = Vq sm{ujt + (f) in [eV] . The spaces 
before and after the middle point are treated as drift 
space. The transverse velocities are modified to keep the 
transverse momentum unchanged. 



4-. Transverse quadrupole resonator 

The TR field is handled in two parts. Longitudinal 
components at the ends are modeled as a longitudinal 
kick. The transverse component inside the TR are mod- 
eled similarly to the magnetic quadrupole. The equa- 
tion of motion is derived starting with dp/dt = Eg using 
eq. (H)) and taking into account that the 7 is changing 
during the transfer, we obtain 

x"{t)-f{t) + x'{t)^{t) - ^ cos(0 + ujt)x{t) = , 
niQ 

y"{th{t) + y'mt) + ^ cos(0 + ojt)y{t) = , 
mo 



z"{th{t) + z'mt)^o , 



where 



7(i)3 {x'{t)x"{t) + v'{t)y"{t) + z'{t)z"{t)) 



The same equation was derived from the relativistic La- 
grangian as a cross check. We ignored the magnetic field 
inside the rods in this model. There is some weak mag- 
netic field between the rods, but if we ignore both this 
weak magnetic field and the small longitudinal variation 
of the transverse electric field our model stays consistent. 

At the entry and the exit of the TR there are longitu- 
dinal components of the EM field which must be taken 
into account. These end fields are modeled as longitudi- 
nal kicks. The procedure is very similar to the one used 
for the accelerating cavity with some difference. The en- 
ergy change is not only a function of the particle's phase 
respect to the synchronous particle, but also a function 
of the transverse position at the entry and at the exit. 
The energy change is given by eq.(I3]). 



B. Tracking results 

As we have already seen, the force introduced by the 
TR can not produce damping. When the end fields of 
the TR were included, the damping of the betatron os- 
cillations and the increase of the synchrotron oscillation 
amplitude described in the first attempt of this paper 
has disappeared. The transverse and longitudinal emit- 
tances stayed constant. There are however special con- 
ditions when the emittance can be transferred between 
the planes. When one of the transverse tune and the 
synchrotron tune satisfy the condition 



nq±^ = mqs 



(16) 



(14) 



synchro-betatron resonance occurs. Where n,m are in- 
tegers and q is the fractional part of the tune. Usu- 
ally the the coupling is introduced by some device, like 
a dipole mode cavity in a dispersive region. We found 
that a transverse-longitudinal coupling is always present 
without the need for any special device. Even when we 
simulated the simplest FODO cell consisting of a 

focusing and a defocusing magnetic quadrupole and an 
accelerating cavity, the synchro-betatron resonance was 
present, if the tunes were adjusted to satisfy the reso- 
nance condition in eq. (|16p . The path length of a particle 
always depends to some extent on its betatron motion. 

A particle with small betatron amplitude has a shorter 
path length than a particle with a large betatron ampli- 
tude. We made simulations with machines containing the 
TR, then the TR was replaced by a magnetic quadrupole 
giving the same focusing strength. There was no notice- 
able difference between the behavior of the two system. If 
the TR introduce some coupling between the transverse 
and longitudinal planes it is much smaller than the intrin- 
sic coupling, which is present even without it. What is 
the role of the TR after all? The emittance tends to flow 
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FIG. 3: Single particle phase-space plots showing synchro- 
betatron resonance between the horizontal and longitudinal 
planes. The ellipses are fitted to points at different stages 
during the emittance transfer. 



FIG. 4: Initial distributions of 9621 points in phase-space and 
after 1300 turns. 



from the plane with bigger emittance to the the smaUer. 
This process has to be stopped at a certain moment to be 
useful. In a dedicated ring is proposed for the emit- 
tance transfer tuned on synchro-betatron resonance. The 
beam is injected into that ring, then it makes a few hun- 
dred turns while the emittance transfer is taking place, 
then it is ejected. The TR can be useful for emittance 
transfer in existing rings. Its role is to modify the trans- 
verse tunes rapidly to satisfy the condition in eg. pop . 



then move the transverse tunes away from the synchro- 
betatron resonance. With normal magnetic quadrupoles 
this can not be done, they are too slow. The voltage of 
the TR can be changed fast enough to do that. There 
is an other possibility to introduce resonance conditions 
quickly by changing the voltage of the accelerating cav- 
ity. This way the longitudinal tune is changed. This 
second method has a serious disadvantage. The longi- 
tudinal tune is usually a small value, in our simulations 
it was 0.024. This means the transverse tune has to be 
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very close to an integer value. With a transverse tune so 
close to integer resonance the beam would suffer losses 
during a longer period of time. Using the TR the trans- 
verse tune can be set to this dangerous value only for the 
duration of the emittance transfer, in our simulation it 
was around 1300 turns. 

Fig-® shows singe particle simulation results in a ring 
at 1.1 GeV/c. The tunes are qh = 7.0244, = 7.281, qi = 
0.024. These tunes were obtained with four TR in the 
ring , one in each quadrant giving 11.8 MV/m electric 
field. These tunes correspond to the case when m = n — 
1 in eq. p^ . The ellipses are fitted to the phase space 
points at different times. As it can be seen on Fig.® the 
emittance start to flow from the bigger to the smaller 
one. It is possible to drive the beam from the smaller 
emittance to the bigger one when n ^ m in ea. (|16p . but 
when n or m bigger than 1 the resonance range is very 
narrow. The transverse tunes have a small amplitude 
dependence, this leads to smearing of the phase space as 
it is noticed in Q. 

When the TRs were off, the horizontal tune was far 
from resonance, qh = 7.1. 

Fig.® shows bunch tracking results. Emittances were 
calculated before and after tracking with the method 
given by Wolski [l^. His method can be used even 
if there are correlations between the planes . Here is a 
short summary. The beam S matrix is defined as = 
{{X,~X;)iXj~X])), where = (x, x', y', z, ^) 
and a matrix S defined as 
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Then the eigenvalues of SS" are , ±iey , ize^. 

The brackets () mean averaging over all particles. The 
result of the emittance calculations are in TABLEU) As it 
can be seen, the 6D rms emittance is slightly bigger after 
tracking. The RMS emittance can grow even the system 
is symplectic. We think it is due to the nonlinear nature 
of the coupling between the transverse and longitudinal 
planes. The product of the three eigenvalues, which is 
the 6D RMS omittance is approximately conserved. 
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FIG. 5: The evolution of emittances as a function of turns. 



About 1300 turns later the horizontal and longitudinal 
planes have the maximum amount of correlation, at this 
point the two emittances come closest to each other. An 
other 1300 turns later the emittances are nearly back to 
the initial conditions, then the cycle starts again. 



V. SUMMARY AND OUTLOOK 

We have revised the idea presented in the first attempt 
of this paper and re-cited the mechanism preventing us 
from getting non-Hamiltonian effects. Maxwell's equa- 
tions imply eq.®, then it follows that the two damp- 
ing terms in the equations of synchrotron motion are 
equal and compensate each other. The force acting on 
the beam is conservative. We found there is an inherent 
coupling between the transverse and longitudinal planes. 
This coupling is present even when there is no any spe- 
cial device in the ring to provoke it. We proposed to use 
the TR to put the tunes of a ring quickly to synchro- 
betatron resonance for the duration of the emittance 
transfer, then to move the tunes away from resonance. 
With our method there is no need for a separate ring. 
It can have applications, where transferring the bigger 
emittance from one plane towards the plane with smaller 
emittance has an advantage. 



VI. 
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Eh 




e; 




initial 


39.3 


9.87 


5.02 


1.95E-9 


final 


16.4 


9.86 


13.8 


2.24E-9 



TABLE I: Initial Rms emittances and final emittances after 
1300 turns . Emittances (la) are in [mm mrad]. 

Fig.® shows the evolution of the emittances calcu- 
lated from the E matrix as mentioned before. The ini- 
tial distributions are Gaussian and without correlations. 
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